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SUMMARY 

An engineering analysis and computer code (AERSEP) for predicting NASA 
Shuttle Orbiter - HO Tank longitudinal aerodynamic characteristics during 
abort separation has been developed. Computed results are applicable at 
Mach numbers above 2 for angle-of-attaek between + 10 degrees. No practical 
restrictions on orbiter-tank relative positioning are indicated for tank 
under-orbiter configurations. Input data requirements and computer running 

times are minimal facilitating program use for parametric studies, test 

I 

planning, and trajectory analysis. In a majority of cases AERSEP Orbiter- 
Tank interference predictions are as accurate as state-of-the-art estimates 
for interference-free or isolated-vehicle configurations. AERSEP isolated- 
orbiter predictions also show excellent correlation with data. 
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INTRODUCTION 


Wind tunnel tests to define Shuttle Orbit er- Tank abort staging character- 
istics are particularly time consuming and expensive. In addition to the many 

! 

test variables the problem is compounded by the extensive data required to 
define the aerodynamic characteristics during the early most critical phase 
of tank separation. Given these circumstances it is imperative to maximize 
the information obtained during test by early identification of major problem 
areas. Recognizing the need NASA MSC awarded Grumman the present contract to 

I 

develop an jengineering capability for predicting Shuttle Orbiter-Tank 
longitudinal characteristics during abort separation. The primary flight 
envelope was to cover Mach No. >2 and angle- of- attack between 
+ 10 degrees. Secondary objectives were to develop a corresponding capability 
for large angle-of-attack, > 25°, and to assess the feasibility of a mid 
angle-of-attack analysis, 10° < < 25°. 

Both primary and secondary contract objectives have been successfully 
accomplished. This report documents the engineering analysis for the 
Aer odynamics of Sep aration, AERSEP, computer code satisfying the primary 
contract objective. The code prepared by Mr. A. Vachris of the Grumman 
Aerodynamics Department is functionally structured to permit program 
generalization as future needs may dictate. AERSEP is compatible with the 
NASA JSC UNTVAC 1108 Executive 8 System. Input data requirements and running 
time are minimal facilitating its use for parametric studies, test planning, 
and trajectory analysis. 

Technical requirements dictated a non-linear analysis which would 
reproduce the predominant orbiter-tank aerodynamic interference phenomenon 
during staging. The selected approach combines principal elements of 
slender body theory, linear theory, and accepted hypersonic approximations. 



The data comparisons in Section 10 of this report indicate AERSEP Orbiter-Tank 
interference predictions are consistent in accuracy with state-of-the-art 
computer predictions limited to interference-free or isolated vehicle situations 
AERSEP isolated orbiter predictions also show excellent correlation with data. 

Results to date indicate that the present AERSEP code can be generalized 
to cover the mid and high angle- of-attack range. Accomplishing this would 
require added program logic to phase in the appropriate engineering analysis 
best suited to each angle-of-attack range. This phasing can become complicated 
since different criteria would apply to forebodies, afterbodies, and wings. 
Developing L satisfactory mid o< interference analysis may have to be 
predicated on the availability of isolated orbiter wing-body and body-alone 
test data. The high 0< case could be patterned after the analysis of 
Mr. R. Hendrikson of Grumman already incorporated into a computer code 
delivered to NASA MSC under the present contract. For the near term it is 
recommended that AERSEP be modified to accept the glove geometry character- 
istizing the current Shuttle Orbiter double-delta design.Until such time the 
present code can be employed in a manner that will approximate the net glove 


effect. 
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Section 1: General Discussion 
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Sections 2 through 8 of the present report document the engineering 
analysis j for the AERSEP code. Section 9 summarizes input data requirements 
and offers practical recommendations on configuration modeling. This section 
provides a brief overview emphasizing the key assumptions, approximations , and 
practical expediencies adopted to realize a working code within contract 
resource limitations. 

The major objective is to develop an engineering capability to predict 

i 

Shuttle Orbjiter/HO Tank longitudinal aerodynamic characteristics during 
staging/abort- separation. The primary flight envelope covers M > 2 and a 
between + 10 degrees. The crux of the problem is to evaluate the orbiter- 
tank mutual interference since isolated configuration data are generally available. 
It is apparent that aerodynamic estimates for the mated vehicle will inherently 
be less accurate than isolated vehicle predictions since both the flow field 
and the vehicle loading must now be calculated. Experience suggests that 
NASA MSC 040-A mated vehicle interference estimates within A ^ A S .010 to 
.015 (based on reference wing area and MAC) are consistent with the accuracy 
level of such programs as the GAC-HAPP or MacDonnell-Douglas HARB for isolated 
vehicles. Larger discrepancies would indicate the potential improvement 
available through further refinement of either the analysis or computer code. 

A review of AERSEP correlation with NASA test data indicate that the majority 
of AERSEP estimates satisfy the stated criteria. 

The analysis divides into two essential steps; determining ^ j 

the combined flow field of the orbiter-tank combination and evaluating the 
surface pressure loading on both vehicles. A rigorous solution for either is 
beyond the state-of-the-art. To break the deadlock AERSEP employs an iterative 



I 


approach. The first approbation assumes that the incident flow on the 
orbiter is the same as the isolated tank flow field. 

I 

The analytic basis for constructing the isolated tank and orbiter flow 

i 

I 

fields is covered in Sections 2, 3, and 4. Perhaps the most notable feature 

I * 

here is the selection of a Mach dependent non-linear body flow field analysis 
which accounts for fuselage fineness ratio , camber, and cross-sectional 
shape effects. The predicted flow fields approach the correct slender body 
limit for high fineness ratio configurations at low M, low o(, ^^^c^/O^axid 
the appropriate hypersonic limit when non-linear effects predominate. 

The next step is to evaluate the "Primary Loading" on the orbiter as 
situated in the isolated tank flow field. Section 5 describing the fuselage 
primary load evaluation may appear simple, however, success depends critically 
on the technique employed. The analysis in AERSEP will approximate linear 
theory results at low Mach number (for slender configurations) and approaches 
the appropriate hypersonic limit at high Mach number (M>5). Both linear and 
non-linear contributions are accounted for at intermediate Mach No., and 

fineness ratio. 

Primary loads on the orbiter wing are obtained by determining the trace 
of the tank shock/wing intersection and integrating the non-uniform upwash 
loading over the affected wing area. Orbiter fuselage -on -wing upwash is 
allowed for by including a interference factor (Section 6). Primary 

loads for the HO Tank are calculated in the same manner as for the orbiter. 

The last phase of the aerodynamic analysis, "reflected waves," is covered 
in Sections 7 and 8. It was noted earlier that the isolated tank flow field 
serves as the first approximation to the combined flow field as seen by the 
orbiter. As a result the primary loads on the orbiter do not include the 



/ 

I 

I. 

effect of orbiter disturbances reflecting off the tank and back onto the 
orbiter. , AERSEP accounts for this by constructing an improved approximation. 

For the orbiter the incident flow is now taken to be the sum of the isolated 

I 

i 

tank flow field (first approximation) and a perturbation flow due to the primary 
loads on the tank. The strength of the outgoing distrubances, i.e., tank 
reflections, is determined from slender body theory and their outward propagation 
calculated along the characteristic structure of the isolated tank flow field. 

The subsequent (improved) loads on the orbiter are then calculated in much the 
same manner! as the primary loads. 

At this point the aerodynamic loading on the tank and orbiter in proximity 
have been evaluated. To obtain the corresponding interference characteristics 
the isolated orbiter and tank loading must also be calculated. Isolated orbiter 
force and moment characteristics are handled as a special case by the AERSEP 
primary loading subroutines. 

Comparison of AERSEP estimates and wind tunnel data for the isolated tank 

shows a consistent tendency to underestimate normal force and pitching moment. 

This is attributed to tank forebody-on-afterbody lift carryover not accounted 

for in the current AERSEP code. Although the net effect of this carryover on 

orbiter-tank interference estimates should be less (compensating effects) 

their inclusion in AERSEP is one potential area for refinement. Wing Wist 

and camber are also not accounted for in AERSEP since their net effect on 

orbiter-tank interference is small compared to root incidence and angle-of- 

attack. The only consequence is that C estimates for the isolated orbiter 

m 0 

are invalid. AERSEP isolated longitudinal stability estimates, however, 
show excellent agreement with NASA MSC 040-A data at M = I .96 and 4.0 using 
= 1.0 (see Figure 10.3). 
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AERSEP aerodynamic and geometric input data requirements are treated in 

Section 9 * It may "be noted that only a minimum of "set up time" is required. 

| 

Computer running time is low, approximately 9 CPU seconds for 6 vertical 

F 

separations at a given M,^< , and incidence on the IBM 370/165. 



Section 2: Isolated Flowfield - HO Tank 

The Shuttle Orbiter wing-body flowfield in AERSEP is represented by an 

i 

overlay of the isolated fuselage and exposed wing flows. Since the general 

j 

approach is simply an extension of that for the HO Tank it is expedient to 

i 

begin with the isolated tank analysis. 

The tank consists of a blunted nose cone, Sketch 2-A, followed by a 
constant diameter cylindrical section with a short boat-tail. 



c fW^T'c/J : /Vas£ Cos/£ 


Ignoring tip blunting the nose flow is conical throughout the region ABC, 
AB'C' determined by the closing characteristics (Mach waves BC and B'C') 
defining the upstream influence of the shoulder. All flow variables in this 
region can be obtained from the exact inviscid flow tabulations in 
Reference 1 and 2 for unyawed/ yawed pointed cones. For obvious reasons these 
numerical data had to be condensed before inclusion into the AERSEP code, 
specifically into four correlation plots. 


? 



The ' windward centerline surface pressure is given by the expression 


^ " ' 71 . T 


4 - A 
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Eq. 2:1 


where 


Cp = cone surface pressure coefficient ato<f = 0 


lsCj? - increment in surface pressure due to angle of attack, <=< 

2 . 

The zero angle-of-attack term is obtained from Figure 2-1 where 0 C 

( 0 = cone half angle) is correlated with the supersonic-hypersonic similarity 
parameter.. /£, 77// 77 A/ © c > . 

The numerical data: in this figure are from Ref. 1. Note that for large 
Tan © c the exact cone solutions approach the Newtonian result = 

7? (Tbs ^' . For moderate or low ^>Tan © c , however, the 9" at c>7 - 0 is 
substantially greater than given by Newtonian. AERSEP calculations are based 
on the faired curve shown. 

The second term in Eq. 2:1 accounts for angle-of-attack effects on the 
nose cone surface pressure. Despite the fact that Newtonian theory does not 
consistently give good results for the surface pressure atoC = 0 it 
adequately predicts the change in surface pressure with angle-of-attack. 

This is amply demonstrated in Figure 2-1. Accordingly, 


A c g<* = ^ ^ 'ST/ 
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Eq. 2:2 




Similarly, for the leeward centerline 



Eq. 2:3 




The corresponding nose cone shock wave strength is determined from the 
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equat ion 
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Eq. 2:4 


where 

dyO 


static pressure coefficient immediately behind wave ato<^ = 0 


The wave to surface pressure ratio at cx^ = 0 appearing in the first term is 
given by the correlation in Figure 2-3. This plot also demonstrates that most 
of the cone surface pressure rise at low ^/3/d/V is ih e result of 

isentropic compression through the cqne shock layer. Conversely, shock 
compression is the dominant mechanism at large • Angle-of-attack 

effects on are given by the second term in Eq. 2:4 where the ratio 

7 ^, /&b(J is given in Figure 2-4. The selected AERSEP 
fairing favors are the most significant Shuttle cases, i.e., . ~ 25°, 

oC ^ - io°. 

Once Cyp is known the shock wave angle, , (relative to the freestream) 

can be calculated from the equations: 
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where 

S 

S* 


component of freestream Mach No. normal to shock wave 

flow deflection through shock 

Mach No. immediately behind shock wave 


/Z 



At I this point all essential flow properties at the shock and surface are 
known with the exception of the surface Mach No . , /^ . The latter is related 
to /%/ by the isentropic relation 



'-Flow-properties’within ‘the shock layer are evaluated assuming a linear 

variation with distance between the cone surface and shock. 

| 

A special provision is included in AERSEP to preclude the predicted pressure 


from falling below the base pressure level ^ event 

lower pressures are predicted the code over-rides the prediction with a 


^ //A/ statement. This requirement is justified by available high speed 


orbiter pressure data. 


The HO tank flow development downstream of the nose cone is characterized 
by an expansion fan emanating from the shoulder as shown in Sketch 2-B. In 
AERSEP the overall fan is subdivided into small expansion fan increments which 
describe the progressive surface flew turning from ^ (relative to the tank 
FRL) on the nose to 0° on the cylindrical tank section. In two dimensional 
flow the required Mach No. variation with expansive flow turning angle can be 
determined from the well known Prandtl-Meyer relation. This equation, however, 
is not generally applicable in three dimensional flow and gives completely 
erroneous results near the shock wave at low supersonic Mach Nos. In AERSEP 
the problem is resolved by visualizing the expansion fan as a succession of 
conical flows generated by a stacked set of cone frustrums as in Sketch 2-C 


below. 




\ \ 

Assum ing , as in Shock-Expansion Theory, that .reflected wave effects . 

are negligible it follows that the flow properties along the outgoing 
characteristic from "Frustrum 2" are identical to those for the conical flow 
field generated by a pointed cone of half angle • This means that the 

surface and shock wave properties for each frustrum can be determined in 
exactly the same manner as for the nose cone. For the HO tank in Sketch 2-B 



the cone-cylinder shoulder is treated as a limiting case of the geometry in 
Sketch 2-C, i.e., B g , B 3 coalesce at the comer. 

A substantially different analysis must be employed to determine the 
flow field generated by the cylindrical section of the tank where^^/V^. "O 
AERSEP' s description in this region is based on yawed cylinder reasoning. 


The prescribed centerline surface variation with 'K. is given by the 

following equations. 


Lower Surface 
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Eq. 2:11 
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Upper Surface 



A little patience will verify that the foregoing equations give the 
variation shown in Sketch 2-D(a) and 2-D(b) for the upper and lower centerlines 
respectively. The resulting . varies as shown in Sketch 2-D(c). 

Interpreting and as average pressures acting over the full width 

of the tank cylinder , it follows that the normal force per unit length of tank 

cylinder is 


. 4 / 




For C 5<L ' > ‘ ' > i.e., for supersonic crossflow, Newtonian theory and test 





data 


predict =0.67 for a circular cross-section. For 

z 





.5 to . 65 . 


i, e. j subsonic crossflow, the data of Ref. 4 indicates 
For practical purposes one can assume the crossflow is constant for all 

M and cx^. The only significant flowfield distinction between the subsonic 
and supersonic crossflow case is that leeward base drag is predominant in the 
former while windward pressure drag (shock wave drag) dominates in latter. 


This distinction is recognized in AERSEP by the surface formulation given 
in Eqs. 11 through 18. 

The corresponding shockwave properties for the tank-cylinder are given 
by Eqs. 2:20 to 2:28 below. 
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surface Mach Ho . , />/j , for each of the above can be determined from 
the isentropic relation 



Eq. 2:24 
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Upper surface Mach No . , , is again calculated from /^/ and Eq. 2:24 
Discussion of the tank boat-tail and base expansion flowfields will be 
deferred to the next section where arbitrary bodies (non-axisyrametric and 
cambered) are considered. 
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Section 3*. Isolated Flowfield - General Body Shapes 

An engineering method for calculating the isolated flow field surrounding 
the Shuttle Orbiter fuselage is given below. The general approach is a 

I 

direct extension of that for the isolated tank and includes the effect of 
non-circular cross-sections and fuselage camber. The required geometry data 
are indicated in Sketch 3-A. 

The orbiter fuselage "true" inboard profile in Sketch 3-A is taken 
directly from the configuration lines drawing. AERSEP approximates the true 
profile by J cutting the fuselage into a series of "plugs" with straight line 
connects between controlling fuselage stations, e.g., FS 155? 250? ^50? 600, 
800, etc. The following geometrical input data are required at the selected 

fuselage control stations. 

- Z to upper surface of true inboard profile 

= Z to lower surface of true inboard profile 

^ — z to upper surface of "effective" inboard profile (see Sketch 3~B) 

n: = ^ to lower surface of "effective" inboard profile (see Sketch 3-B) 

iV = Fuselage width 

These data are used by AERSEP to calculate the following derived geometric 
information. 

= Cross-sectional area at designated fuselage control station 


sf = . tX Y'Ze u - 


Eq. 3:1 


= Z distance to centroid of cross-sectional area at designated 


fuselage control stations 
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Eq. 3=2 



& Cfj ~ Equivalent cone angle for n'th fuselage plug 
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Eq. 3:3 


<£<- = Equivalent cone incidence for n'th plug 
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Eq. 3:^ 


Eq. 3:5 


d? = Upper surface incidence of n'th plug 
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Eq. 3:6 


The total fuselage flow is constructed by overlaying the individual flows 
generated by each fuselage plug. Each plug in turn is considered to be a 
yawed equivalent cone (at angle-of-attack and incidence) which has been 
distorted to conform to the actual fuselage cross-sectional. The expression 
for surface pressure is accordingly made up of three terms; 


do 
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Eq. 3:7 
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where 


Co- 


i Surface pressure coefficient on "equivalent" cone at oL - 0 


A L increment to surface pressure coefficient due to equivalent 

l cone angle -of- attack plus incidence = cxf 


A C 0 - Increment to surface pressure coefficient due to distorting the 
XP, 

^ "equivalent" circular cone to the actual fuselage cross-sectional 
shape. 

These are evaluated as follows: 
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from Figure 2-1 
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Eq. 3:8 




Eq. 3:9 


(upper, surface ) 
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The rationale and justification underlying the and _A 

contributions was treated earlier (Eq. 2:l). Justification for the A 
term as given in Eq. 3'10 and 3 ‘-12 paralles that for the A 
contribution, Eq. 3:8 and 3:9- In each case it is assumed that the effect of 
local parcel angle-of-attaek changes due to the angle -of- attack, * 

incidence, ; and distortion to a non-circular cross-section can be 

evaluated using Newtonian Theory. Figure 2-2 substantiates this assumption 
for effec ^ s 8X1(1 inference for C * effects. Figure 3-1 and 

3-2 provide additional justification and demonstrate the suitability of 
Eqs. 3:7 thru 3:13 as a general method for evaluating the orb iter fuselage 
surface pressure distribution. 

A similar procedure is used to determine the shock wave properties 
corresponding to a fuselage "plug". The pressure coefficient immediately 
behind the shock is given by 
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Eq. 3:1^ 
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where 

(Z-/;? _ = Surface pressure coefficient behind shock wave of "equivalent 

i 

cone at = 0 

I 

= Increment to shock wave pressure due to equivalent cone angle-of- 
attack plus incidence, =<. 


-A. .Co = Increment to shock wave pressure due to distorting the "equivalent 
^ circular cone to the actual fuselage cross-sectional shape. 

These terms are evaluated as follows: 
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Eq. 3=15 


(lower surface) 
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Eq. 3=16 




(upper surface) 
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Eq. 3:17 
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The correlation curve in Figure 3~3 for the ratio of shockwave to surface 
shape change was obtained by multiplying the correlation in Figure 2-4 by the 

(^/^/V Z ’ ZisS- J/./ 




factor 
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I 

This formulation agrees with slender body theory at low M which predicts the 

j 

radial perturbation velocity due to cross-sectional shape change will decay 

l/r times more rapidly than the radial velocity due toO<^ • As M & the 

j 

above factor approaches unity which is the appropriate limit for hypersonic 

flow With C/s> and known the remaining flow properties (3 ^ , > 

<S ^ 

•/^and can be determined from Eqs. 2:5 to 2:10 given earlier. 

Expansion fans originating from the juncture of two successive fuselage 
plugs (at a fuselage control station) are built up in the same manner as the 
shoulder expansion fan for the HO tank as described in that section. The only 
new aspect in the present case is that each expansion fan has an L c and 
associated with it which are determined by linear interpolation between 
adjacent plug values. 

The preceding analysis applys to all fuselage plugs for which 
The general fuselage analysis when 0 — 0.15 is a direct 

generalization of Eqs. 2:11 to 2:18, Eqs 2:20 to 2:26, and Eqs. 2:27 to 2.30 
for the isolated tank. Including incidence and boat-tail effects these 
equations now read. 
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The corresponding shock wave properties are given below 
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The AERSEP formulation for the tank and orbiter base flowfield at zero 

angle of attack is based on the simple model shown in Sketch 3-C. The surface 

Mach No., M 0 , and flow direction, , immediately upstream of the tank 

S/ Js. 


base are Mg^= M ^ 


- V 

and . (Rel to FRL) = 0 if, .as done here, .we 
neglect the relatively small boat-tail. Based on simple linear theory the 
wake streamline direction is given by 
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where 






Eq. 3:42 


The local pressure, , and flow direction, <£§2^, anywhere behind 


the Mach wave off the tank shoulder are related to the corresponding wake 


streamline values by 
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The base flow analysis for the orbiter lower surface at OC~" -O is 
virtually identical, i.e., 




and 


Assuming the base pressure is constant for a < 10° it follows that 

(Rel to Free-Stream) are inariant with <*. The above power-off base flow 
representation is invalid beyond, say, one diameter aft of the base at M = 2. 

AERSEP calculations are inapplicable when far -wake pressures can exert 

i . . 

significant interference effects. Any doubt can be resolved by a simple 
order-of-magnitude analysis accounting for Mach wave sweep, area affected, and 
interference pressure level determined from Eqs. '3‘. ^6 and 3:^7- 
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Section 4: Isolated Wing- Body Flow Field 


For the intended flight envelope, M — 2, — + 10°, the wing leading 

edge will he supersonic. In effect a Mach wave drawn from the exposed root 
chord leading edge at o<C = 0 lies aft of the leading edge and the upper and 
lower surface flows do not interact except for a small region near the tip. 
With increasing angle of attack the wing flow, viewed in a plane perpendicular 
to the leading edge, will appear as shown in Sketch 4-A. For a thin planar 
wing an oblique shock forms at the underside of the leading edge to redirect 
the flow parallel to the surface. As increases the deflection angle £ 
will eventually equal the shock detachment angle £, 1 / 4 / corresponding to the 
normal Mach No. , PL. Any further increase in ©< results in leading edge 
shock detachment .and significant windward- leeward cross flow, i.e., three 
dimensional flow. Below this critical c< an essentially two dimensional flow 
field description is adequate. Sketch 4-A gives the variation in o^-for-shock- 
detachment with for various . 

A prof ile view of the AERSEP lower wing surface flow is shown in Sketch 
4-B. The flow deflection due to the wing is ' 


■ <§>iy ~~ ^ E l* 

where is the exposed root chord incidence and the body upwash 

factor. Wing thickness, camber, and twist are neglected. When ** 0 the. 

; streamwise shock wave angle (Eel. to freestream) can be determined by 

iterative solution of the equation. 
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Mach Wo. behind the wing shock, , is calculated from 


= /C 


Eq. 4:3 




Eq. 4:4 


The corresponding surface pressure is given by 
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Eq. 4:5 


At negative orbiter angle-of-attack the bottom of the wing can become 
a leeward surface, i.e., ^ 0. For the small negative angles of interest 

the pressure can be calculated from the linearized expression 
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Eq. 4:6 


(also require ^ alwa y s ) 
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Eq. 4:7 


Since the Shuttle wing trailing edge is coincident with the orbiter base, 
the flow at the trailing edge will re-expand to the base pressure level 






z, 

'CO 


. For 4, 2 0 the Prandtl-Meyer relations at the 


trailing edge are written as 



Eq. 4 : 8 


■= /Or " <% j 


where <==^r , jf? and >/* are the flow direction (relative to FRL, 

measured positive up) and pressure on the wing lower surface and a point in 
the expansion fan respectively. Mach No. M, at a point in the fan is given 
'hy 



To further simplify AERSEP logic the shock wave curvature introduced hy the 
trailing edge expansion waves is neglected and the shock wave location 
approximated hy a straight line. 

When 0 the trailing edge expansion fan relations are 



An additional feature included in AERSEP to partially account for three 
dimensional orhiter wing-hody interference effects is illustrated in Sketch 
4-C. Spillover of the wing lower surface pressure onto the fuselage is 
confined to a region aft of the Mach waves from the root chord leading 



edge. For supersonic leading edges linear theory predicts the average wing 
carry-over pressure, , across the fuselage width, W, will vary from zero 

at the root leading edge, Station "A", to the exposed wing pressure level at 
Station "C" located t wo fuselage Mach wave crossings aft of "A". Theory also 
indicates the longitudinal variation is approximately linear as indicated in 
the lower part of Sketch 4-C. AERSEP approximates this distribution over the 
orbiter fuselage with a "step" loading at Station "B", one Mach-wave- crossing 
aft of the exposed root leading edge. 


Sketch 4-D gives a planview of the orbiter wing-body overlaid on the 
tank. The hatched region indicates the extent of tank plan area blanketed by 
wing pressures. For simplicity this irregular area is converted to an 
equivalent rectangle extending from B* to the orbiter base, i.e., 



fZ}-" t ’/) / ) 



The net orbiter lower surface wing-body flow field is thus given by the 
isolated fuselage flow field everywhere forward of the wing shock from Station B* 
and by the two dimensional wing flow downstream. A representative AERSEP 
calculation for the NASA 040-A orbiter is plotted to scale in Sketch 4-E. 
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Section 5: Primary Interference Loads on Bodies 


Perhaps the most critical step in the entire AERSEP analysis is the 
evaluation of body loads in a non-uniform interference flow field. "Primary" 
refers to the aerodynamic loading on the orbiter when placed in the isolated 
tank flow field. "Reflected" interferences due to orbiter shockwave 
reflections off the tank back onto the orbiter are accounted for separately. 

The geometry for the primary load evaluation on the orbiter fuselage is 
indicated in Sketch 5-A. The tank and upper half of its flow field are shown 
underneath the "effective" outline of the orbiter. As discussed in Section 3 
the "effective" outline measures from the FRL to the average fuselage height 
as seen in cross-section. Local flow properties evaluated along the "effective' 
orbiter outline are more representative for load evaluation than those along 
the true geometric profile. 

The average primary interference pressure across the fuselage width of 
the orbiter lower surface at some station X is given by 








Eq. 5:1 


where 






//= is 




f 

/t~ 








~ Ratio, cross-flow drag coefficient of orb iter fuselage lower 
surface to that of a flat plate 

~ A cross-flow drag factor that accounts for the quasi -ax isymmetric 
nature of the tank interference flow field; a function of 
W = fuselage width at station X 

= Angle-of -attack between orbiter "effective" lower surface at 
station X and the local flow direction in the tank interference 
flow field 

= Local static pressure and Mach No. in the tank flow field along 
orbiter "effective" lower surface 

4 . 


To evaluate Eq. 5:1 is calculated from 


ex'* « ^ 


Eq. 5:2 




where 0 is the flow direction in the tank flow field relative to the tank 
f 

FRL. Both and are measured relative to their respective FRL's. 

# 

The orbiter effective lower surface incidence, Ca , is given by 

= / ^Zj// ~ / Eq. 5:3 

vC- / •yy - . -- J 

where "n" denotes "at the n'th" fuselage control station. Appropriate values 
for and the decay constant "C" (which determines K) can be read-from 

Figure 5-1 and 5-2 for each fuselage plug. Average upper "effective" surface 
loading is determined in exactly the same manner and the resulting force and 
moment determined by numerical integration. 



Isolated orbiter fuselage force and moment are calculated in a similar 
manner. The appropriate equations for the orbiter lower surface are: 


/ co 
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Eq. 5:^ 
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-Eq. 5:5 
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Eq. 5 : 6 


Equation 5:5 for is nothing more than the familiar Newtonian expression 

for surface pressure coefficient multiplied by a cross-flow drag ratio to 

account for the cross-sectional shape. 

In light of the opening remarks to this section it might seem unlikely 
that Eq. 5:1 would be adequate for load evaluation over the Mach number range 
2 to 10. In point of fact, however, the integrated force and moment predicted 
by this equation are identical to slender body theory when is small and 

Newtonian -Cross flow results when nonlinear effects are significant. These 

theoretical anchor points are vital and assure reasonable AERSEP predictions 
over the intended Mach spectrum. ; 
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Section 6: Primary Interference Loads, on Wing-Body Configurations 


The total lift on the Shuttle Orbiter in free air can be expressed as 


-S 


^3 




Eq. 6:1 


where 

,<S - 

c 




Lift of isolated body at OC 

Lift curve slope of exposed wing panels when joined together to 
form an isolated wing * 

= Interference factor which accounts for the lift increment on the 
exposed wing due to fuselage induced upwash at ' :> <C 
~ Interference factor which accounts for lift carry over from the 
exposed wing at<^ onto the fuselage 
^*33 ~ Interference factors due to wing incidence; similar to 

For configurations with supersonic leading edges and no fuselage aft of 
the wing trailing edge Eq. 6:1 can be simplified to 
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Eq. 6:2 


^2^] depends on the body width to wing span ratio, w/b, as shown in 
Figure 6-1 from Ref. 5» 

Equation 6:2 can be adapted to the present case for evaluating the orbiter 
loading in the interference flow field generated by the tank. Primary 
interference on the fuselage, corresponding to > can calculated by 



the method obtained in the previous section. The additional loading on the 
orb iter due to wing and wing -body contributions is determined from the 
loading integral 
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where 


A c, 
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t/L. 





Local dynamic pressure 

Differential pressure coefficient across wing surface referred to 
the local dynamic pressure 

Total angle- of -attack between wing surface and the local flow 

direction in the tank flow field 

Orbiter wing incidence (relative to orbiter FKL) 


The integrand in Eq_« 6j 3 is evaluated assuming that the tank interference 
flow field is approximately axisymmetric (about the wind axis) in the region 
of the wing. Sketches 6-A and 6-B illustrate the procedure for determining 
the off- centerline local flow properties in the wing Mean Reference Plane., 
.MRP. Consider the transverse plane A-A which intersects the centerline shock 
envelope at "S" in Sketch 6-A. A circle' swung through "S" as in Sketch 6-B 
intersects the wing at span station "S*" along W-W, the trace of A-A in 



the MRP. 


Successive calculation of "S*" for a range of fuselage stations 
defines the shock trace on the exposed wing. Flow properties in each meridian 
plane, 0 = constant, are assumed ,to be the same as in the plane of 
symmetry. The total angle- of- attack, , between the local flow and the 

wing at some point "P*" is calculated from 
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Eq. 6:4 


In performing the load integration, Eq. 6:3, the exposed wing is divided 
into regions internal and external to the shockwave trace. The tank flow field 

d. 

defines the local flow within the trace while free stream coqgitions prevail 
outside. The wing contribution to is given by 
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Section 7: Wave Reflections from Bodies 


It will be recalled from Sections 5 and 6 that the "Primary Interference" 
loads determined thus far do not account for shockwave reflections. A 
representative situation is illustrated in Sketch. 7-A. In general the 
strength of the reflected flow field will be proportional to the primary 
interference loading and, for bodies, will decay with distance from the 
reflecting surface. Reflected intererences will be subordinate to the primary 
loading and slender body theory should be adequate for its evaluation. 

For simplicity consider the special case of orbiter wave reflections off 
the HO tank at = 0. Let* the primary interference loading on the 

tank over the length AX. at station X/ be shown in Sketch 7-B and zero elsewhere. 
According to slender body theory the potential function for the perturbation 
flow field associated with this incremental load is 


Q 7Z.r) = - a-X) Cxsft 

- r 


Eq. 7:1 


where X , K ' are the cylindrical coordinates of an arbitrary field point. 
Equation 7:1 describes the induced flow due to a lineal dipole singularity 
distribution of strength stretched out along the tank X axis (FRL). 

The local strength factor can be evaluated by equating the dipole 

normal force to the tank primary interference loading over the interv al . X/ 
to -X/ . Thus 
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I'he normal force due to the primary interference loading can be written as 



Hence, from Eq. 7:2 and 7:3 



Eq. 7:3 




7:4 


Substituting back into Eq. 7:1 for the potential it follows that 




The left; hand side of Eq. 7:6 may be revitten as follows 
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Eq. 7:7 


where fa is the average reflected <fa> , i.e., Cy? over the interval X./ 

' A 

to . The final expression for is 
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On the tank surface in the plane of symmetry, j ^ ~~ & j _ 




Eq. 7:9 


For 

i.e. 


fa fa the reflected wave strength varies inversely with /"* , 
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Simple slender body theory (above) implies that the flow in the crossflow 
plane through X/ is completely determined by the body cross-sectional shape 
'dnd loading at Xy • In actual fact reflections off the tank surface will be 
swept aft along the local Mach wave as shown in Sketch 7-B. The AERSEP code 
carries this one step further and propagates the tank reflections along the 
characteristic structure (Mach wave net) of the isolated tank flow field. 

Once the reflected wave is known the remaining reflected flow 

properties are determined from the relations 



Eq. 7:11 



Eq. 7:12 



Eq. 7:13 


It should be noted that <^0^ > the average flow direction defined above, 
is in a sense fictitious. The reflected disturbance from the tank at X/ 
influences the entire flow field aft of the Mach cone from Xy and is not con- 
fined to the AX wide ray indicated in Sketch 7-B. The net effect of this 
assumption results in only minor error in the calculated linear reflected 



are 


loading since the integrated "Sf and the centroid of its distribution 
well predicted. Nonlinear loads, which depend on local 6^, , would be 
unreliable, however, and are ignored in the AE.RSEP reflected load calculation. 
Order of magnitude analysis " indicates these non-linear loads will be small. 

Tank shockwave' and flow field reflections off the orbit-er fuselage are 
handled in the same way as reflections off the tank. The appropriate equations 
for the orbiter, corresponding to Eq. 7: 9 and 7:10 for the tank, are 



where '2^- , 2"^ , 2 are the 2 . coordinates of the orbiter lower surface, 

orbiter camber line, at Z/ and a field point on the same orbiter character- 
istic (Mach ray) as 2/ , the surface point. 

Reflected loads are determined using the same AERSEP computer routines as 
for the primary interference load calculation where, now, the reflected flow 
field is substituted for the isolated flow field. The only significant change 
is a modification to Eq. 6:4 which, for reflected loads on the wing, now _ 

•reads . 
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Eq. 7: l6 



1 '‘J.X . 

The added cos factor is due to the cos (j6 dependence of the reflected 


dipole flow field. 




Section 8: Wave Reflections from Wing-Bodies 


The reflected flow field from the orbiter wing lower surface is defined 
in much the same way as the isolated wing flow field in Section 4. The 
primary interference load or normal force on the exposed wing panel may be 
written as 


if 


Eq. 8:1 


■ './////§ jex/w&p ~ 

In Eq. 8:1 and ' are the differential.^. across the wing 

( lower minus upper) in the isolated tank flow field and the 

freestream respectively. These terms can be identified with the loading 
integral, Eq. 6:3. Approximately one half of the primary interference loading 
in Eq. 8:1 is due to the primary loading on the lower surface. Thus, the 
refected wave strength from the wing onto the tank can be approximated by 
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Eq. 8:2 


The integral in Eq. 8:2 extends over the inboard exposed wing root chord. 

The reflected flow field geometry is taken to be identical to the isolated 
orbiter wing-body flow field. Orbiter fuselage reflections define the flow 
back to the wing shock as in Sketch 4-B. The wing reflected flow prevails 
aft of the shock emanating from point "B" in Sketch 4-C. 
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Section 9: Summary and AERSEP Input Data Requirements 

AERSEP input data requirements are modest as may be seen by the following 
MSA M3Ci040-A data-set. Sketch 9-A portrays the 040-A orb iter "True", 
"Approximate", and "Effective" profiles (see Section 3 for definitions) used 
to generate the AERSEP predictions appearing in Sections 10.2 to 10.7. Vertical 
ordinate data and fuselage width at each control station are indicated in 
Sketch 9-B. Effective crossflow drag ratios, , and decay constants, C, 

for each fuselage plug are also noted. Both parameters are determined by the 
local cross-sectional shape using Figures 9-1 and 9-2 except in nose regions 

I. ’ . 

where. Cp/g. =1.0. In addition AERSEP’ s fuselage representation is 
stereotyped as constant-section from the exposed root chord leading edge aft. 

Input data for the wing is limited to a specification of the leading 

edge location (Fuselage Station) at the exposed root and tip, the reference 

* 

wing semi -span, and wing root incidence, . The present code assumes a 

simple trapazoidal planform with an unswept trailing edge terminating at the 
fuselage base. Where some (modest) fuselage extension beyond the wing and/or 
trailing edge sweep is encountered ignore the extended fuselage length and 
substitute an equivalent exposed wing panel with the same root chord, area, and 
5O/ 0 MGC location. The effects of orbiter fuselage upwash on the wing is 
accounted for by a factor (see Figure 6-l) which, within reason, may 

alternately be selected to match isolated orbiter predicted with data. 

s HO Tank data appear in Sketch 9-C which parallels the orbiter fuselage 
presentation. Both orbiter and tank moment reference center locations are 
specified in their individual body axis system. 



AERSEP nomenclature and sign convention for the relative orb iter -tank 
positioning is indicated in Sketch The primary data are 2r , <=> ^, , 

and . Tank angle of attack, Ck v 7"" , is a derived quantity. More detailed 
information concerning the -AERSEP code is given in the users manual. Reference 7 
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j Section 10: Data Comparisons - General Results 

| 

AERSEP predictions for representative configurations are compared with wind 

/ • -- 

tunnel data in Figures 10.1 through 10.7* In all figures the experimental data 
are denoted by open symbols and AERSEP predictions by heavy crosses. All test 
data have been plotted to as large a scale as warrented by reading accuracy from 
the original source. 

AERSEP results and data for parallel staged tanks at M = 2,5 and 10 are 
compared in Figures 10.1-1 through 10.1-9. The zero stagger side by side 
arrangement results in high levels of interference. For the closest spacing at 

M = 2 the note shock from each impinges on the forebody of the other resulting 

! ' 

in extreme surface pressures and undoubtedly some local efessfe? reflected shock 
detachment. Despite this the interference moment is reasonably well predicted 
at all o<^ . The normal force correlation is less satisfactory. 

Pitching moment and normal force predictions at M = 5 show good agreement 
with data barring the normal force ato{ =8? At M = 10 an extreme sensitivity 
to small variations in separation distance is noted and the AERSEP results are 
slightly displaced from the available data. It is clear, however, that the 
heavy nose-down moment at the close spacing must reverse to a strong nose-up 
as the tanks separate and the impinging shock "walks" back. The high M results 
imply that the analysis is sensitive to the precise selection of receiver 
"effective height" at which average incident flow properties for determining 
Tank "2" loads are evaluated. The need to select an "average" point can be 
avoided by employing a three dimensional receiver surface subdivided into small 
panels for local loads. evaluation. This alternative was considered when form- 
ulating the present code but was ruled out by schedule limitations. 



/ 


AERSEP calculations and isolated OUO-A Tank data at M - 1.96 and. ^.0 are 
compared in Figures 10.2-1, 10.2-2. The predicted results are independent of 
Mach No. , ’underestimate the lift at low angle -of-attack, and underestimate the 
pitching moment at M/ = l.'96.“ ' In Common with" contemporary engineering programs 
AERSEP ignores fuselage forebody-on-afterbody lift carry-over effects. The 
carry-over lift aerodynamic center is near the shoulder at low Mach No. and 
moves aft at high ' Mach No. which ■accounts for the pitching moment discrepancy 
at M = 1.96.' Nonetheless the agreement with data indicates that AERSEP s 
loading analysis is adequate for the present purposes. 

Isolated 040-A orbiter characteristics and AERSEP predictions are com- 
pared in Figure 10.3-land 10.3-2. AERSEP does not account for wing twist or 
camber effects, hence, the Cj^ prediction should be ignored. This simplification 
has virtually no effect on predicted interference levels and does not compromise 
AERSEP 's utility for this purpose in practical abort -separation cases. Lift 
curve slope and longitudinal stability show excellent correlation with data at 
M = I .96 and M = 4.0. 

Representative NASA 040-A Orbiter/Tank wind tunnel data (Ref. 6 ) and AERSEP 
interference predictions are compared in Figures 10.4 through 10.7. No simple 
generalization of these data is possible except to note that in a majority of 
cases the AERSEP results are within Ac^, A C M - .015 of the test data, i.e., 
within the accuracy level of contemporary arbitrary body computer programs limited 

to interference-free configurations. _ - - - • 

Earlier in this section it was noted that the present code uses a minimum of 

information to describe the Orbiter/Tank geometry. Instead AERSEP employs cross 
flow drag coefficients, radial flow factors, a body upwash factor, and "effective 
fuselage upper/ lower heights to determine average values which summarize the 
principal configuration effects. Preliminary 040-A AERSEP correlations indicated 

Ss ~A 



that an alternative technique for calculating average flow direction in the 

I 

tank (Orhiter) afterbody flowfield gave* better results than the method pre- 
scribed in Section 3 and 5* Instead of using the flow direction in the 
vertical plane of symmetry (Eq. 3:29-3:39) the average direction is now 
related to the local C . Spacifically where the generator /&A/ @ c “ ' & 
the average local flow direction is given by: 


Lower Surface 


Upper Surface 


yo /TTer&Jr/T&TT/J — 





With the advantage of hindsight the explanation is almost obvious. At angle 
of attack the subsonic crossflow velocity distribution around, say, the tank 
will vary betweeen twice the freestream value at the maximum breadth to zero 
(rel. to FRL) in the vertical plane of symmetry. Moving away from the tank 
surface there is a strong tendency for the average crossflow to return to the 

freestream value, ^(rel to freestream) « c> \/^c 5 S ° that in effect at 
small + C*(_ the adjacent orbiter fuselage is not shielded by the belly mounted 


tank. 

With regard to the selection of other AERSEP parameters the following 
should be noted. Extensive experience with delta wing shuttle orbiter 
configurations has consistently shown that selecting - 1 gives the best 

isolated orbiter predictions. The same proves true for predicting interference 
effects. Slightly better correlation was also observed when the AERSEP wing 
shock was defined to originate from the exposed root chord leading edge instead 
of the B* location prescribed in Section 4. This choice partially compensates 
for wing thickness effects not included in the AERSEP. analysis. 



/ 

I 
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In reviewing various 040-A AERSEP calculations and laying out the computed 

/ 

flowfield structure it was noted that at no time did the HO Tank lower surface 
reflections impinge onto the orbiter. The code was modified accordingly so 
that only orbiter-lower-surf^ce/tank-upper-surface primary loads were used to 
determine the reflected wave strength. 

All AERSEP results in this section were calculated using the procedural 
modifications noted above. In summary the data correlations suggest that 
extending the present AERSEP code to incorporate a three dimensional model of 
the orbiter-tank geometry would result in consistently better estimates 
particularly at MS: 4. As it is, by judicious averaging, the present AERSEP 
code prediction accuracy (with interference ) approaches that expected for 
interference-free situations. ■ ' ■ 


*? 



CONCLUSION 


M engineering analysis and AERSEP computer code for predicting NASA 

Shuttle Orbit er - HO Tank longitudinal aerodynamic characteristics during 

I 

abort separation has been developed. Computed results are applicable at 
Mach number above 2 for angle -of -attack between + 10 degrees'. ' No practical 
restrictions on orbiter-tank relative positioning are indicated for tank- 
under-orbiter configurations. AERSEP is also compatible with the NASA JSC 
UNIVAC 1108 Executive System. Input data requirements and running time 
are minimal] facilitating its use for parametric studies, test planning, 
and trajectory analysis. In a majority of cases AERSEP Orbiter-Tank 
predictions with interference are comparable in accuracy with state-of-the- 
art estimates for interference-free configurations. 

Further generalization and/or refinement of the present AERSEP form- 
ulation to expand the applicable flight envelope is feasible. Extending 
AERSEP to incorporate a detailed three dimensional orbiter-tank geometry 
model would result in improved estimates throughout the present flight 
envelope M~2, o( + 10°. 
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